library(irr)

# matriz com todos os arquivos
codif_todos <- read.csv2("codif.csv")
codif_todos_m <- as.matrix(codif_todos[,c(6:8)])

## ka
ka_todos <- kripp.alpha(t(codif_todos_m),
                         method = "nominal")

## percentual identico nos 3
perc_id3_todos <- sum(rowSums(codif_todos_m)==3|
        rowSums(codif_todos_m)==0)/nrow(codif_todos_m)

## percentual igual aos pares
perc_id2_todos <- ((sum(rowSums(codif_todos_m[,c(1,3)])==2|
        rowSums(codif_todos_m[,c(1,3)])==0)/nrow(codif_todos_m))+
(sum(rowSums(codif_todos_m[,c(1,2)])==2|
        rowSums(codif_todos_m[,c(1,2)])==0)/nrow(codif_todos_m)))/2



# matriz sem repets
codif_sr <- subset(codif_todos, dupl==0)
codif_sr_m <- as.matrix(codif_sr[,c(6:8)])

## ka
ka_sr <- kripp.alpha(t(codif_sr_m),
                         method = "nominal")

## percentual identico nos 3
perc_id3_sr <- sum(rowSums(codif_sr_m)==3|
        rowSums(codif_sr_m)==0)/nrow(codif_sr_m)

## percentual igual aos pares
perc_id2_sr <- ((sum(rowSums(codif_sr_m[,c(1,3)])==2|
        rowSums(codif_sr_m[,c(1,3)])==0)/nrow(codif_sr_m))+
    (sum(rowSums(codif_sr_m[,c(1,2)])==2|
           rowSums(codif_sr_m[,c(1,2)])==0)/nrow(codif_sr_m)))/2

# tabela
round(
  data.frame(
  todos=c(ka_todos$value, perc_id3_todos, perc_id2_todos),
  sem_repet=c(ka_sr$value, perc_id3_sr, perc_id2_sr),
  row.names = c("krip_alfa", "perc_ident_todos", "perc_ident_pares")
), 2)

